Ancestry dynamics and trait selection in a designer cat breed.

The Bengal cat breed was developed from intercrosses between the Asian leopard cat, Prionailurus bengalensis, and the domestic cat, Felis catus, with a last common ancestor approximately 6 million years ago. Predicted to derive ∼94% of their genome from domestic cats, regions of the leopard cat genome are thought to account for the unique pelage traits and ornate color patterns of the Bengal breed, which are similar to those of ocelots and jaguars. We explore ancestry distribution and selection signatures in the Bengal breed by using reduced representation and whole-genome sequencing from 947 cats. The mean proportion of leopard cat DNA in the Bengal breed is 3.48%, lower than predicted from breed history, and is broadly distributed, covering 93% of the Bengal genome. Overall, leopard cat introgressions do not show strong signatures of selection across the Bengal breed. However, two popular color traits in Bengal cats, charcoal and pheomelanin intensity, are explained by selection of leopard cat genes whose expression is reduced in a domestic cat background, consistent with genetic incompatibility resulting from hybridization. We characterize several selective sweeps in the Bengal genome that harbor candidate genes for pelage and color pattern and that are associated with domestic, rather than leopard, cat haplotypes. We identify the molecular and phenotypic basis of one selective sweep as reduced expression of the Fgfr2 gene, which underlies glitter, a trait desired by breeders that affects hair texture and light reflectivity.

NDUFS7 variant in dogs with Leigh syndrome and its functional validation in a Drosophila melanogaster model.

Two Jack-Russell Terrier × Chihuahua mixed-breed littermates with Leigh syndrome were investigated. The dogs presented with progressive ataxia, dystonia, and increased lactate levels. Brain MRI showed characteristic bilateral symmetrical T2 hyperintense lesions, histologically representing encephalomalacia. Muscle histopathology revealed accumulation of mitochondria. Whole genome sequencing identified a missense variant in a gene associated with human Leigh syndrome, NDUFS7:c.535G > A or p.(Val179Met). The genotypes at the variant co-segregated with the phenotype in the investigated litter as expected for a monogenic autosomal recessive mode of inheritance. We investigated the functional consequences of the missense variant in a Drosophila melanogaster model by expressing recombinant wildtype or mutant canine NDUFS7 in a ubiquitous knockdown model of the fly ortholog ND-20. Neither of the investigated overexpression lines completely rescued the lethality upon knockdown of the endogenous ND-20. However, a partial rescue was found upon overexpression of wildtype NDUFS7, where pupal lethality was moved to later developmental stages, which was not seen upon canine mutant overexpression, thus providing additional evidence for the pathogenicity of the identified variant. Our results show the potential of the fruit fly as a model for canine disease allele validation and establish NDUFS7:p.(Val179Met) as causative variant for the investigated canine Leigh syndrome.

PCYT2 deficiency in Saarlooswolfdogs with progressive retinal, central, and peripheral neurodegeneration.

We investigated a syndromic disease comprising blindness and neurodegeneration in 11 Saarlooswolfdogs. Clinical signs involved early adult onset retinal degeneration and adult-onset neurological deficits including gait abnormalities, hind limb weakness, tremors, ataxia, cognitive decline and behavioral changes such as aggression towards the owner. Histopathology in one affected dog demonstrated cataract, retinal degeneration, central and peripheral axonal degeneration, and severe astroglial hypertrophy and hyperplasia in the central nervous system. Pedigrees indicated autosomal recessive inheritance. We mapped the suspected genetic defect to a 15 Mb critical interval by combined linkage and autozygosity analysis. Whole genome sequencing revealed a private homozygous missense variant, PCYT2:c.4A>G, predicted to change the second amino acid of the encoded ethanolamine-phosphate cytidylyltransferase 2, XP_038402224.1:(p.Ile2Val). Genotyping of additional Saarlooswolfdogs confirmed the homozygous genotype in all eleven affected dogs and demonstrated an allele frequency of 9.9% in the population. This experiment also identified three additional homozygous mutant young dogs without overt clinical signs. Subsequent examination of one of these dogs revealed early-stage progressive retinal atrophy (PRA) and expansion of subarachnoid CSF spaces in MRI. Dogs homozygous for the pathogenic variant showed ether lipid accumulation, confirming a functional PCYT2 deficiency. The clinical and metabolic phenotype in affected dogs shows some parallels with human patients, in whom PCYT2 variants lead to a rare form of spastic paraplegia or axonal motor and sensory polyneuropathy. Our results demonstrate that PCYT2:c.4A>G in dogs cause PCYT2 deficiency. This canine model with histopathologically documented retinal, central, and peripheral neurodegeneration further deepens the knowledge of PCYT2 deficiency.

Single-cell profiling of bronchoalveolar cells reveals a Th17 signature in neutrophilic severe equine asthma.

Severe equine asthma (SEA) is a complex respiratory condition characterized by chronic airway inflammation. It shares many clinical and pathological features with human neutrophilic asthma, making it a valuable model for studying this condition. However, the immune mechanisms driving SEA have remained elusive. Although SEA has been primarily associated with a Th2 response, there have also been reports of Th1, Th17, or mixed-mediated responses. To uncover the elusive immune mechanisms driving SEA, we performed single-cell mRNA sequencing (scRNA-seq) on cryopreserved bronchoalveolar cells from 11 Warmblood horses, 5 controls and 6 with SEA. We identified six major cell types, including B cells, T cells, monocytes-macrophages, dendritic cells, neutrophils, and mast cells. All cell types exhibited significant heterogeneity, with previously identified and novel cell subtypes. Notably, we observed monocyte-lymphocyte complexes and detected a robust Th17 signature in SEA, with CXCL13 upregulation in intermediate monocytes. Asthmatic horses exhibited expansion of the B-cell population, Th17 polarization of the T-cell populations, and dysregulation of genes associated with T-cell function. Neutrophils demonstrated enhanced migratory capacity and heightened aptitude for neutrophil extracellular trap formation. These findings provide compelling evidence for a predominant Th17 immune response in neutrophilic SEA, driven by dysregulation of monocyte and T-cell genes. The dysregulated genes identified through scRNA-seq have potential as biomarkers and therapeutic targets for SEA and provide insights into human neutrophilic asthma.

Independent COL17A1 Variants in Cats with Junctional Epidermolysis Bullosa.

Epidermolysis bullosa (EB), characterized by defective adhesion of the epidermis to the dermis, is a heterogeneous disease with many subtypes in human patients and domestic animals. We investigated two unrelated cats with recurring erosions and ulcers on ear pinnae, oral mucosa, and paw pads that were suggestive of EB. Histopathology confirmed the diagnosis of EB in both cats. Case 1 was severe and had to be euthanized at 5 months of age. Case 2 had a milder course and was alive at 11 years of age at the time of writing. Whole genome sequencing of both affected cats revealed independent homozygous variants in COL17A1 encoding the collagen type XVII alpha 1 chain. Loss of function variants in COL17A1 lead to junctional epidermolysis bullosa (JEB) in human patients. The identified splice site variant in case 1, c.3019+1del, was predicted to lead to a complete deficiency in collagen type XVII. Case 2 had a splice region variant, c.769+5G>A. Assessment of the functional impact of this variant on the transcript level demonstrated partial aberrant splicing with residual expression of wildtype transcript. Thus, the molecular analyses provided a plausible explanation of the difference in clinical severity between the two cases and allowed the refinement of the diagnosis in the affected cats to JEB. This study highlights the complexity of EB in animals and contributes to a better understanding of the genotype-phenotype correlation in COL17A1-related JEB.

RALGAPA1 Deletion in Belgian Shepherd Dogs with Cerebellar Ataxia.

Several genetically distinct forms of cerebellar ataxia exist in Belgian shepherd dogs. We investigated a litter in which two puppies developed cerebellar ataxia. The clinical signs stabilized at around six weeks of age, but remained visible into adulthood. Combined linkage and homozygosity mapping delineated a 5.5 Mb critical interval. The comparison of whole-genome sequence data of one affected dog to 929 control genomes revealed a private homozygous ~4.8 kb deletion in the critical interval, Chr8:14,468,376_14,473,136del4761. The deletion comprises exon 35 of the RALGAPA1 gene, XM_038544497.1:c.6080-2893_6944+1003del. It is predicted to introduce a premature stop codon into the transcript, truncating ~23% of the wild-type open reading frame of the encoded Ral GTPase-activating protein catalytic subunit α 1, XP_038400425.1:(p.Val2027Glnfs*7). Genotypes at the deletion showed the expected co-segregation with the phenotype in the family. Genotyping additional ataxic Belgian shepherd dogs revealed three additional homozygous mutant dogs from a single litter, which had been euthanized at five weeks of age due to their severe clinical phenotype. Histopathology revealed cytoplasmic accumulation of granular material within cerebellar Purkinje cells. Genotyping a cohort of almost 900 Belgian shepherd dogs showed the expected genotype-phenotype association and a carrier frequency of 5% in the population. Human patients with loss-of-function variants in RALGAPA1 develop psychomotor disability and early-onset epilepsy. The available clinical and histopathological data, together with current knowledge about RALGAPA1 variants and their functional impact in other species, suggest the RALGAPA1 deletion is the likely causative defect for the observed phenotype in the affected dogs.

SGCD Missense Variant in a Lagotto Romagnolo Dog with Autosomal Recessively Inherited Limb-Girdle Muscular Dystrophy.

An 8-month-old female Lagotto Romagnolo dog was presented for a 1-month history of an initial severe reluctance to move, rapidly progressing to a marked stiff gait and progressive muscular weakness and evolving to tetraparesis, which persuaded the owner to request euthanasia. A primary muscle pathology was supported by necropsy and histopathological findings. Macroscopically, the muscles were moderately atrophic, except for the diaphragm and the neck muscles, which were markedly thickened. Histologically, all the skeletal muscles examined showed atrophy, hypertrophy, necrosis with calcification of the fibers, and mild fibrosis and inflammation. On immunohistochemistry, all three dystrophin domains and sarcoglycan proteins were absent. On Western blot analysis, no band was present for delta sarcoglycan. We sequenced the genome of the affected dog and compared the data to more than 900 control genomes of different dog breeds. Genetic analysis revealed a homozygous private protein-changing variant in the SGCD gene encoding delta- sarcoglycan in the affected dog. The variant was predicted to induce a SGCD:p.(Leu242Pro) change in the protein. In silico tools predicted the change to be deleterious. Other 770 Lagotto Romagnolo dogs were genotyped for the variant and all found to be homozygous wild type. Based on current knowledge of gene function in other mammalian species, including humans, hamsters, and dogs, we propose the SGCD missense variant as the causative variant of the observed form of muscular dystrophy in the index case. The absence of the variant allele in the Lagotto Romagnolo breeding population indicates a rare allele that has appeared recently.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture.

BACKGROUND: The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function. RESULTS: We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection. CONCLUSIONS: We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Case report: Sacral agenesis in two boxer dogs: clinical presentation, diagnostic investigations, and outcome.

Two boxer dogs from the same litter were presented at 3 months of age for urinary and fecal incontinence. Both dogs had an abnormal tail consisting of a small stump, an atonic anal sphincter, and absent perineal reflex and sensation. Neurological evaluation was indicative of a lesion of the cauda equina or sacral spinal cord. Radiology and CT scan of the spine displayed similar findings in the two dogs that were indicative of sacral agenesis. Indeed, they had 6 lumbar vertebrae followed by a lumbosacral transitional vertebra, lacking a complete spinous process, and a hypoplastic vertebra carrying 2 hypoplastic sacral transverse processes as the only remnant of the sacral bone. Caudal vertebrae were absent in one of the dogs. On MRI, one dog had a dural sac occupying the entire spinal canal and ending in a subfascial fat structure. In the other dog, the dural sac finished in an extracanalar, subfascial, well-defined cystic structure, communicating with the subarachnoid space, and consistent with a meningocele. Sacral agenesis-that is the partial or complete absence of the sacral bones-is a neural tube defect occasionally reported in humans with spina bifida occulta. Sacral agenesis has been described in human and veterinary medicine in association with conditions such as caudal regression syndrome, perosomus elumbis, and Currarino syndrome. These neural tube defects are caused by genetic and/or environmental factors. Despite thorough genetic investigation, no candidate variants in genes with known functional impact on bone development or sacral development could be found in the affected dogs. To the best of the authors' knowledge, this is the first report describing similar sacral agenesis in two related boxer dogs.

PKD1 Nonsense Variant in a Lagotto Romagnolo Family with Polycystic Kidney Disease.

A female Lagotto Romagnolo dog with polycystic kidney disease (PKD) and her progeny, including PKD-affected offspring, were studied. All affected dogs appeared clinically inconspicuous, while sonography revealed the presence of renal cysts. The PKD-affected index female was used for breeding and produced two litters with six affected offspring of both sexes and seven unaffected offspring. The pedigrees suggested an autosomal dominant mode of inheritance of the trait. A trio whole genome sequencing analysis of the index female and her unaffected parents identified a de novo heterozygous nonsense variant in the coding region of the PKD1 gene. This variant, NM_001006650.1:c.7195G>T, is predicted to truncate 44% of the open reading frame of the wild-type PKD1 protein, NP_001006651.1:p.(Glu2399*). The finding of a de novo variant in an excellent functional candidate gene strongly suggests that the PKD1 nonsense variant caused the observed phenotype in the affected dogs. Perfect co-segregation of the mutant allele with the PKD phenotype in two litters supports the hypothesized causality. To the best of our knowledge, this is the second description of a PKD1-related canine form of autosomal dominant PKD that may serve as an animal model for similar hepatorenal fibrocystic disorders in humans.

Autosomal recessive hyposegmentation of granulocytes in Australian Shepherd Dogs indicates a role for LMBR1L in myeloid leukocytes.

Pelger-Huët anomaly (PHA) in humans is an autosomal dominant hematological phenotype without major clinical consequences. PHA involves a characteristic hyposegmentation of granulocytes (HG). Human PHA is caused by heterozygous loss of function variants in the LBR gene encoding lamin receptor B. Bi-allelic variants and complete deficiency of LBR cause the much more severe Greenberg skeletal dysplasia which is lethal in utero and characterized by massive skeletal malformation and gross fetal hydrops. HG phenotypes have also been described in domestic animals and homology to human PHA has been claimed in the literature. We studied a litter of Australian Shepherd Dogs with four stillborn puppies in which both parents had an HG phenotype. Linkage analysis excluded LBR as responsible gene for the stillborn puppies. We then investigated the HG phenotype in Australian Shepherd Dogs independently of the prenatal lethality. Genome-wide association mapped the HG locus to chromosome 27 and established an autosomal recessive mode of inheritance. Whole genome sequencing identified a splice site variant in LMBR1L, c.191+1G>A, as most likely causal variant for the HG phenotype. The mutant allele abrogates the expression of the longer X2 isoform but does not affect transcripts encoding the shorter X1 isoform of the LMBR1L protein. The homozygous mutant LMBR1L genotype associated with HG is common in Australian Shepherd Dogs and was found in 39 of 300 genotyped dogs (13%). Our results point to a previously unsuspected function of LMBR1L in the myeloid lineage of leukocytes.

Ancient segmentally duplicated LCORL retrocopies in equids.

LINE-1 is an active transposable element encoding proteins capable of inserting host gene retrocopies, resulting in retro-copy number variants (retroCNVs) between individuals. Here, we performed retroCNV discovery using 86 equids and identified 437 retrocopy insertions. Only 5 retroCNVs were shared between horses and other equids, indicating that the majority of retroCNVs inserted after the species diverged. A large number (17-35 copies) of segmentally duplicated Ligand Dependent Nuclear Receptor Corepressor Like (LCORL) retrocopies were present in all equids but absent from other extant perissodactyls. The majority of LCORL transcripts in horses and donkeys originate from the retrocopies. The initial LCORL retrotransposition occurred 18 million years ago (17-19 95% CI), which is coincident with the increase in body size, reduction in digit number, and changes in dentition that characterized equid evolution. Evolutionary conservation of the LCORL retrocopy segmental amplification in the Equidae family, high expression levels and the ancient timeline for LCORL retrotransposition support a functional role for this structural variant.

A TNR Frameshift Variant in Weimaraner Dogs with an Exercise-Induced Paroxysmal Movement Disorder.

BACKGROUND: Some paroxysmal movement disorders remain without an identified genetic cause. OBJECTIVES: The aim was to identify the causal genetic variant for a paroxysmal dystonia-ataxia syndrome in Weimaraner dogs. METHODS: Clinical and diagnostic investigations were performed. Whole genome sequencing of one affected dog was used to identify private homozygous variants against 921 control genomes. RESULTS: Four Weimaraners were presented for episodes of abnormal gait. Results of examinations and diagnostic investigations were unremarkable. Whole genome sequencing revealed a private frameshift variant in the TNR (tenascin-R) gene in an affected dog, XM_038542431.1:c.831dupC, which is predicted to truncate more than 75% of the open read frame. Genotypes in a cohort of 4 affected and 70 unaffected Weimaraners showed perfect association with the disease phenotype. CONCLUSIONS: We report the association of a TNR variant with a paroxysmal dystonia-ataxia syndrome in Weimaraners. It might be relevant to include sequencing of this gene in diagnosing humans with unexplained paroxysmal movement disorders. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

SOAT1 missense variant in two cats with sebaceous gland dysplasia.

Spontaneously arisen hereditary diseases in domestic animals provide an excellent opportunity to study the physiological functions of the altered genes. We investigated two 4-month-old sibling domestic short haired kittens with dry dark debris around the eyes, nose, and ears, dark crusting on the legs and a thin poor hair coat. Skin biopsies revealed abnormal sebaceous gland morphology with lack of normal sebocyte arrangement and differentiation. Hair follicles had a distorted silhouette, interpreted as a change secondary to the observed sebaceous gland dysplasia. Whole genome sequencing on both affected kittens and 65 genetically diverse feline genomes was performed. Filtering for variants that were present in both kittens but absent from the control genomes revealed a homozygous missense variant in SOAT1, encoding sterol O-acyltransferase 1. The protein is localized in the endoplasmic reticulum and catalyzes the formation of cholesteryl esters, an essential component of sebum and meibum. The identified SOAT1:c.1531G > A variant is predicted to change a highly conserved glycine residue within the last transmembrane domain of SOAT1, p.Gly511Arg. In mice, variants in Soat1 or complete knockout of the gene lead to the "hair interior defect" (hid) or abnormal Meibomian glands, respectively. SOAT1:c.1531G > A represents a plausible candidate variant for the observed sebaceous gland dysplasia in both kittens of this study. The variant was not present in 10 additional cats with a similar clinical and histopathological phenotype suggesting genetic heterogeneity. SOAT1 variants should be considered as potential cause in hereditary sebaceous gland dysplasias of humans and domestic animals.

Methemoglobinemia, Increased Deformability and Reduced Membrane Stability of Red Blood Cells in a Cat with a CYB5R3 Splice Defect.

Methemoglobinemia is an acquired or inherited condition resulting from oxidative stress or dysfunction of the NADH-cytochrome b5 reductase or associated pathways. This study describes the clinical, pathophysiological, and molecular genetic features of a cat with hereditary methemoglobinemia. Whole genome sequencing and mRNA transcript analyses were performed in affected and control cats. Co-oximetry, ektacytometry, Ellman's assay for reduced glutathione concentrations, and CYB5R activity were assessed. A young adult European domestic shorthair cat decompensated at induction of anesthesia and was found to have persistent methemoglobinemia of 39 ± 8% (reference range < 3%) of total hemoglobin which could be reversed upon intravenous methylene blue injection. The erythrocytic CYB5R activity was 20 ± 6% of normal. Genetic analyses revealed a single homozygous base exchange at the beginning of intron 3 of the CYB5R3 gene, c.226+5G>A. Subsequent mRNA studies confirmed a splice defect and demonstrated expression of two mutant CYB5R3 transcripts. Erythrocytic glutathione levels were twice that of controls. Mild microcytosis, echinocytes, and multiple Ca2+-filled vesicles were found in the affected cat. Erythrocytes were unstable at high osmolarities although highly deformable as follows from the changes in elongation index and maximal-tolerated osmolarity. Clinicopathological presentation of this cat was similar to other cats with CYB5R3 deficiency. We found that methemoglobinemia is associated with an increase in red blood cell fragility and deformability, glutathione overload, and morphological alterations typical for stress erythropoiesis.

SDR9C7 missense variant in a Chihuahua with non-epidermolytic ichthyosis.

Ichthyoses represent a heterogeneous group of cornification disorders that are associated with skin barrier defects. We investigated a 9-month-old Chihuahua showing excessive scale formation. Clinical and histopathological examinations revealed non-epidermolytic ichthyosis and a genetic defect was suspected. We therefore sequenced the genome of the affected dog and compared the data with 564 genetically diverse control genomes. Filtering for private variants identified a homozygous missense variant in SDR9C7, c.454C>T or p.(Arg152Trp). SDR9C7 is a known candidate gene for ichthyosis in humans and encodes the short-chain dehydrogenase/reductase family 9C member 7. The enzyme is involved in the production of a functional corneocyte lipid envelope (CLE), a crucial component of the epidermal barrier. Pathogenic variants in SDR9C7 have been described in human patients with autosomal recessive ichthyosis. We assume that the identified missense variant in the affected Chihuahua of this study impairs the normal enzymatic activity of SDR9C7 and thus prevents the formation of a functioning CLE, resulting in a defective skin barrier. To the best of our knowledge, this is the first report of a spontaneous SDR9C7 variant in domestic animals.